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Abstract 

The review of recent results in the s=l/2 quantum spin chains with 1/ sinh^(Kr) exchange is 
presented. Related problems in the theory of classical and quantum Calogero-Sutherland-Moser 
systems with inverse square hyperbolic and elliptic potentials are discussed. The attention is 
paid to finding the explicit form of corresponding Bethe-Ansatz equations and to connection 
with generalized Hubbard chains in one dimension. 

1. Introduction 

The idea of spin exchange interaction of electrons as natural explanation of ferromagnetism 
was first proposed by Heisenberg [1] and soon realized in mathematical form by Dirac [2] . 
But the first appearance of the famous Heisenberg Hamiltonian in solid-state physics 
occured three years later in the book by van Vleck [3] . Now it is of common use and was 
investigated from many points of view by various methods of condensed matter theory. 
In two and higher dimensions, the problem of finding the eigenvalues and eigenvectors 
can be solved only by approximate or numerical methods. In one dimension the exact 
solution was obtained in the seminal paper by Bethe [4] who considered most important 
case of nearest-neighbor exchange described by the Hamiltonian 

H= h{j-k){a,ak-l), (1) 

where {crj} are the usual Pauli matrices acting on the s=l/2 spin located at the site j 
and exchange constants {h} are of extreme short-range form, 

/^(J) = ^('^b1,i + 5b1,iv-i). (2) 

It turned out that the solution comes in the form of linear combinations of plane waves 
chosen as to satisfy certain conditions required by (1,2). 

Starting with this solution known as Bethe Ansatz, the investigation of one-dimensional 
exactly solvable models of interacting objects (spins, classical and quantum particles) has 
given a number of results both of physical and mathematical significance [5] . Bethe found 
his solution empirically; at that time the possibility to solve the quantum-mechanical 
problems was not associated with the existence of underlying symmetry. The role of 
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such symmetries has been recognized much later, with one of the highhghts being the 
Yang-Baxter equation which allows one to find some regular way to finding new examples 
of exactly solvable models [6,7]. In many cases, however, the empirical ways are more 
productive since they use some physical information on their background and are not so 
complicated from mathematical viewpoint. 

This concerns especially to the Calogero-Sutherland-Moser (CSM) models which were 
discovered about thirty years ago. They describe the motion of an arbitrary number of 
classical and quantum nonrelativistic particles interacting via two-body singular potentials 
with the Hamiltonian 

M „2 M 

HcsM ^Y.-i + l{l + l)Y. Vi^i - ^k), (3) 

j=l j<k 

where {p, q} are canonically conjugated momenta and positions of particles, I e R and 
the two-body potentials are of the form: 

V{X) = ^, (4) 

V(x) = p(,i), (6) 

where k. G R+ and p{x) is the double periodic Weierstrass p function determined by its 
two periods cui e R+, UJ2 — i'K / k as 



{x — moji — mo-i)^ (mwi + mo^)" 



(7) 



The solvability of the cigcnproblem for first two potentials (4) has been found indepen- 
dently by Calogcro [8] and Sutherland [9] in quantum case while (5) and (6) have been 
found much later [10-11] for classical particles by constructing extra integrals of motion 
(conserved quantities) via the method of Lax pair. Namely, it turned out that the dy- 
namical equations of motion are equivalent to the (M x M) matrix relation 

f = [i^,M], (8) 

where 

ivjfc = Vj^ik + (1 - 3jk)f{xj - Xk), (9) 

M 

Mjk = (1 - 5jk)g{xj - Xk) - 5jk E ^(^i ~ ^'n)' 
if the functions /, g, V obey the Calogero-Moser functional equation 

f{x)g{y) - f{y)g{x) = f{x + y)[V{y) - V{x)] (10) 
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which impUes g{x) = f'{x), V{x) = —f{x)f{—x). The most general form of the solution 
to (10) has been found by Krichcver [12] in terms of the Weierstrass sigma functions 
which give rise to the potential (6). Note that (4) might be considered as hmits of (5,6) 
as K — > 0, and (6) can be regarded as double periodic form of (5) ((5) under periodic 
boundary conditions). The existence of M functionally independent integrals of motion 
in involution follows from the evident relations d{iiL'^) / dt = 0, 1 < n < M. The fact that 
all these conserved quantities are in involution also follows from functional equation (10) 
but needs some cumbersome calculations which can be extended also to the quantum case 
where the time derivative should be replaced to quantum commutator with Hamiltonian 
[13]. In the review paper [13], one can find a lot of interesting facts about the quantum 
models (3-6) established till 1983. 

The Bethe Ansatz technique and the theory of CSM models developed independently 
till 1988 when Haldane [14] and Shastry [15] proposed a new spin 1/2 model with long- 
range exchange resembling (4), 

which has very simple ground-state function of Jastrow type in the antiferromagnetic 
regime J>0 and many degeneracies in the full spectrum. The complete integrabilty of 
the model and the reason of these degeneracies-the sl{2) Yangian symmetry has been 
understood later- for a comprehensive review see [16] and references therein. The Haldane- 
Shastry model has many nice features, including the interpretation of the excited states as 
ideal "spinon" gas, exact calculation of the partition function in the thermodynamic limit 
and the possibility of exact calculation of various correlations in the antiferromagnetic 
ground state [16]. 

The connection with the Bethe case of nearest-neighbor exchange also came soon: in 
1989, I have found that the Bethe and Haldane-Shastry forms of exchange are in fact the 
limits of more general model in which h{j) is given by the elliptic Weierstrass function in 
complete analogy with (6), 

Kj) = JPn{j), (12) 

where the notation pjv nieans that the real period of the Weierstrass function equals N. 
The absolute value of the second period tt/k is a free parameter of the model [17]. The 
Haldane-Shastry spin chain arises as a limit of k — > 0. When considering the case of 
an infinite lattice (A^ — > oo), one recovers the hyperbolic form of exchange (5) which 
degenerates into the nearest- neighbor exchange if k — oo under proper normalization 
of the coupling constant J: J sinh^(/t) J. Hence (6) might be regarded as (5) under 
periodic boundary conditions (finite lattice). Various properties of hyperbolic and elliptic 
spin chains form the main subject of the present review. 

The analogy between quantum spin chains and CSM models is much deeper than 
simple similarity of spin exchange constants and two-body CSM potentials. It concerns 
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mainly in similarities in the form of wave functions of discrete and continuous cases. 
Namely, already in [17] it has been mentioned that the solution of two-magnon problem 
for the exchange (12) and its degenerated hyperbolic form can be obtained via two-body 
CSM systems with potentials (5,6) at Z = 1; it has also been found soon to be true for 
three- and four-magnon wave functions for hyperbolic exchange [18]. Why does this sim- 
ilarity hold? Till now this is poorely understood, but it is working even for elliptic case 
as it will be shown in Sections 3-4. Another question concerns integrability of the spin 
chains with hyperbolic and elliptic exchange, i.e. the existence of a family of operators 
commuting with the Hamiltonian. In the case of nearest-neighbor exchange, such a family 
can be easily found within the framework of the quantum inverse scattering method [6] . 
However, it is not clear up to now how this method should be used in the hyperbolic 
and (more general) elliptic cases. Instead, in Section 2 the Lax pair and empirical way 
of constructing conserved quantities is exposed. Section 5 contains various results for 
hyperbolic models on inhomogeneous lattices defined as equlibrium positions of the clas- 
sical CSM hyperbolic systems in various external fields. Recent results concerning the 
integrability of the related Hubbard chains with variable range hopping are presented in 
Section 6. The list of still unsolved problems is given in the last Section 7 which contains 
also a short summary and discussion. 

2. LcLX pair and integrability 

I shall consider in this Section a bit more general models with the Hamiltonian 

'^AT = ^ Yl hjkPjk, (13) 
l<j¥'k<N 

where {Pjk} are operators of an arbitrary representation of the permutation group S^. 
The spin chains discussed above fall into this class of models, as it follows from the spin 
representation of the permutation group: 

The exchange constants hjk in (13) are supposed to be translation invariant. The notation 
V'jfc = 'ipiJ — k) will be assumed for any function of the difference of numbers j and k 
in this Section. The problem is: how to select the function h so as to get a model with 
integrals of motion commuting with the Hamiltonian (13)? The answer has been done in 
[17]: one can try to construct for the model the quantum Lax pair analogous to (9,10) 
with N X N matrices: 

N 

Ljk = (1 — Sjk)fjkPjk, Mjk = (1 — 5jk)gjkPjk — Sjk X] hjsPjs- 
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The quantum Lax relation \H, L] = [L, M] is equivalent to functional Calogero-Moser 
equation for /, g, h: 

fpqQqr Qpqfqr fpr{hqr flpq 

) (14) 

supplemented by the periodicity condition 

\q — ^p,q+Ni (1^) 

where h'{x) is an odd function of its argument, h^^ = fqpQpq — fpqQqp- The most gen- 
eral solution to (14) has been given in [12] as the combination of the Weierstrass sigma 
functions. There is the normalization of / and h which allows one to write the relations 

h{x) = f{x)f{-x\ gix) = ^^, h\x) ^^^^^ 



dx dx 
The solution given in [12] looks as 

0'{^ -|— Q/ J 

fi^) = / X / ^ exp(-a;C(a)), h{x) = p{a) - p{x). (16) 
a[x)a[a) 

where 

C'W = -p(x). 

and a is the spectral parameter which does not introduce any new in exchange dynamics. 
The periodicity condition (15) means that all Weierstrass functions in (16) are defined on 
the torus Tjv = C/NZ + i^Z, n e R+ is the free parameter of the model. It is easy to 
see that the exchange (16) reduces in the hmit k — > to the Haldane-Shastry model and 
the limit of infinite lattice size {N — > oo) corresponds to the hyperbolic variable range 
form of exchange. And finally, in the limit k — > oo just nearest-neighbor exchange (2) is 
reproduced as it was already mentioned in the preceding Section. 

However, the problem is not classical one and the existence of the Lax representation 
does not guarantee the existence of the integrals of motion as invariants of the L matrix. 
In fact, just for the problem under consideration the operators trL" do not commute 
with the Hamiltonian. Nevertheless, already in [17] the way of constructing integrals of 
motion on the base of / function of Lax pair was proposed. Namely, it was found that 
the operator 

J{a)= fjkfkifijPjkPki (17) 

commutes with Til Moreover, the dependence of the right-hand side on the spectral 
parameter implies that there are two functionally independent operators bilinear in {P} 
commuting with 7i, 

^1 - E (c(j -k)+ cik -i)+ ai - j))Pj,PM, 
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p'U - k) + p'{k -l) + p'il - j)]P,kPki. 

Very long but straightforward calculations show that Ji^2 mutually commute. 

It turns out [32] that the construction (17) can be generalized for more complicated 
operators with higher degrees of {P}- The basic idea is to use the operators of cyclic per- 
mutations Psi...si = Ps-,s2Ps2S3---Psi_isi and functions F^^..^, = fs^sifsisa-- fsi^jsjsisi which 
are invariant under the action of elements of a group of cyclic permutations of subindices 
(1, ../). If one denotes as $(si, .., si) the functions which are completely symmetric in their 
arguments, and J2ceCi Bsi..si as the sum over all cyclic permutations of the subindices of 
BsL.si, the following properties of the above objects are useful: 

(A) The functions h{x) and h'{x) obey the relation 

ceCs 

(B) The sum Fj^^^^^ = EceQ Fsi..sisi+i does not depend on Si+i. 

(C) The sum Esi^..s,+i ^(si, - ^s,+isi)-Psi..s,+i vanishes for any sym- 
metric function 

(D) The sum 

Sl^..Si 

has a representation in the form sj]^i{^) + sj'^i{^), where 

— ^ y: hs„..,si_„p)j:KsAps,..s^., 

and 

ifl>3, 

The main statement concerning the integrals of motion for the Hamiltonian (13) can be 
proved without the use of the specific form (16) of the solution to the Calogcro-Moscr 
equation. It can be formulated as follows: Let (3 < m < A^) be the linear combinations 
of the operators of cyclic permutations in ordered sequences of N symbols, 

Im=E:^^ E ^^'\s,,..,Sm-2l)F,,..s^_,,Ps,...s^_,,. (18) 

'=U Sl^..Srn-2! 

Then they will give the integrals of motion as it follows from 

Proposition 2.1. The operators commute with Tijv given by (13) if the functions 
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are determined by the recurrence relation 

l<j<k<N;j,k^si,..Sm-2l 



or, equivalently, are given by sums over 21 indices 

I 

.., sm-2i) = {ny E ^o-fe} n Kk.^ 

l<ja<ka<N;{j,k}j^Sl,..Sm-2l <^=^ 



(20) 



where X^jk} equals 1 if the product Yl''a-^p{ja — jf3){ka ~ — kp) differs from zero and 

vanishes otherwise. 

The rigorous proof of the statements A-D can be found in [32]. Here I give only sketch 
of the proof of Proposition 2.1. It is based on the calculation of the commutator 



Sl^..Sn 



(21) 



where $ is symmetric in its variables. With the use of invariance of Fsj^„s„ and -Psi..s„ 
under cyclic changes of summation variables it is easy to show that this commutator can 
be written as 



Jn = n 



4^^ + 4^^ + E(i- 



n — 1 n 
~^ ^2 



- Q ) S.,[n/2] ) 4'} 



where 



~ ^("Sl, .., Sn)-fsi..S„(^SnSn+l ^SlS„+l)-fsi..S„+l) 

Sl^..S„ 

'^nj' '^(•Sl) '5n)-^Sl..Sn(^Si/Sn ~ ^SlSI.+l)-fsl..SI/-fsl.+l..S„■ 

Sl7^■•Sn 

The third term can be transformed with the use of functional equation (14) and cyclic 
symmetry of Psi...s^ and Ps„_,_i..s„ to the form 



~l~-^Si/+i..s„(V'si/+isi..Si, 'psi..s„+i) Psi..s„Psi, 



(22) 



where (fsiS2-Sl+l — fsiS2fs2S3--fsiSl + l9si+lSl- 

Now it is easy to see that the term in the first brackets in (22) disappears due to 
statement (B) and the term in the second brackets vanishes due to the relation 



Fsi..Sl{hsiSl + l ^S;Si + l) ^Sl 



S2..Si + l y^Sj+lSl.-Sp 
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which allows to transform this term to the expression which vanishes upon symmetrization 
in all cyclic changes of Si, ..Siy. Hence the operator (21) contains only cyclic permutations 
of rank (n + 1) and (n — 1). This fact leads to the idea of recurrence construction of the 
operators (18) which would commute with the Hamiltonian. It happens if the functions 
{$} obey the recurrence relation 

Yl ihpsm-2i-i - hps,)^^^\si, ..,Sm-2l-l,p) 

= $('+^)(si, S2, .., Sm-21-2) - $^'+'n«m-2/-l, S2, ■., Sm-2/-2), = 1, 

which can be solved in the form (19) or (20). 

The dependence of (18-20) on the spectral parameter a via the relations (16) allows one 
to conclude that there are several integrals of motion at each m. Namely, the analysis of 
this dependence shows that the operators (18) can be written in the form 

m-2 

where Pm commutes with all operators of elementary transpositions, 7^,^ are linear com- 
binations of {Psi..Sm-2i\ which do not depend on a, and w^_^(q;) are linearly independent 
elliptic functions of the spectral parameter. Nothing is known for the mutual commu- 
tativity of these operators except the explicit result for m = 3 mentioned above. Still 
there is no connection with Yang-Baxter theory, i.e. the corresponding i?-matrix and 
L-operators are unknown. There is an excellent paper by K.Hasegawa [33] which states 
that i?-matrix for spinless elliptic quantum Calogero-Moser systems is Belavin's one, also 
of elliptic type. However, it is not clear how to extend Hasegawa's method to the spin 
case so as to reproduce the rich variety of the operators (18). 



3. The infinite chain 



On the infinite line, the model is defined by the Hamiltonian 

-{BjBk - l)/2. 



H 



■-E- 



K 



2 sinh K(j - k) 



(23) 



where j, /c G Z. At these conditions, only ferromagnetic case J > is well defined. The 
spectrum to be found consists of excitations over ferromagnetic ground state |0 > with 
all spins up which has zero energy. The energy of one spin wave is just given by Fourier 
transform of the exchange in (23), 



e{p) = J 



"2^^(2^^ + 2 



P 



TT 



Ci 



in \ 
2^) 



Ci 



ip_\ 
2kJ 



n 2 




(24) 
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where Weierstass functions pi, Ci are defined on the torus Ti = C/(Z + ^Z), i.e. pi has 
the periods (l,a; = m/n). 

3.1. Two-magnon scattering. The two-magnon problem for the model (23) is 
already nontrivial. One has to solve the difference equation for two-magnon wave function 
ip{ni,n2) which is defined by the relation 

1^ >= V'(^l,^2)SniSn2|0 

where the operator {s^^} reverses spin at the site Ua and {ip > is an eigenvector of the 
Hamiltonian (23). The solution is based on the formula [17] 



^ smh^[a{k + z)] 



cothK(/c + I + z) 



k=—oo 



(yi{z - rp) 



coth(«;/) exp 



pz 



TT 



Ci(a;/2) 



x{pi{z) - pi(rp) + 2 



. / X 2rp /w 



+ 



K, 



-(1 — exp(— ipZ)) 



(25) 



\2j sinh(2K0 

x(Ci(^ + rp) - Ci{z) + 2Ci{rp) - Ci{2rp))}, 

where Vp = —up/ Air and / G Z. 

The proof of (25) is based on the quasiperiodicity of the sum on its left-hand side 
and the structure of its only singularity at the point 2; = on a torus Ti obtained by 
factorization of a complex plane on the lattice of periods (l,cc;). The structure of (25) 
allows one to show that the two-magnon wave function is given by the formula 

gi(pini+p2n2) sinhk(ni - na) + 7] + e^(Pi"2+P2ni) sinhk(ni - ^2) - 7] ,^ , 
^{ni, 712) = —T—7Z —^ ' (26) 



the corresponding energy is 



sinh K{ni — 712) 



where e{pi) are given by (24) and the phase 7 is connected with pseudomomenta pi^2 by 
the relation 



coth 7 = — 
' 2k 



IP2 



2k 



ipi 



2k 



Pi -P2 

TT 



(27) 



This gives, in the limit of k — > 00 (a; — > 0), just the expression for the Bethe phase [4], 
and the additivity of magnon energies takes place. The equation (27) can be rewritten in 
the form 

f{pi) - /(P2) 



where 



coth 7 



fip) - ^Ci 

TT 



2k 



m 
2^ 



ip 
2k 



(28) 



It admits also the representation 



f{p) = iKCoi^- kY^ 

^ n=l 



ip 



coth [ — -\- Kn] -\- coth — — ku 



ip 
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If Pi, 2 arc real, the wave function (26) describes scattering of magnons. The relatively 
simple form of (27) allows one to investigate the bound states of two magnons in detail 
[34]. Namely, in these states the wave function must vanish as |ni — 712] — > oo. It means 
that pi and p2 should be complex with P — Pi + P2 real. The simplest possibility is given 
by the choice 

P P 

Pi = ^ + P2 = — - iq, 

where q is real, and one can always choose g > for convenience. Then vanishing of 
ip{ni,n2) as |ni — —>■ 00 is equivalent to the condition 

( N /(Pi) - /(P2) 1 /on^ 
coth7(pi,p2) = 1^ = 1. (29) 

The structure of the function (28) is crucial for the analysis. It is easy to see that it is 
odd and double quasiperiodic, 

/(p)--/(-p), /(p + 27r) = /(p), /(p + 2z/^) = /(p) + 2«. (30) 

Note that one can always choose q < k due to (30). The equation (29) can be rewritten 
in more detailed form 



FA,) = 1 - ^ 



TT \2k/ \4a 2k/ \4a 2k 



= 0. (31) 



At fixed real P and g, the function (31) is real. Moreover, the relations (30) imply the 
following properties of Fp{q), 

Fp(0) = 1, Fp{q) = -Fp{2k - q), Fp{q) = -Fp{-q) + 2. 

One can immediately see that Fp{k) = but this zero is unphysical: the wave function 
in this point vanishes identically. The physical solution, if exists, must lie in the interval 
< q < K. Such a nontrivial zero exists if the derivative of Fp{q) is positive at q — k, 



fiP l\ fiP 1 



> 0. (32) 



This inequality indeed takes place for the values of P within the interval < P < Per, 
< Per < TT [34]. There should be at least one bound state specified by (31). At P > Per, 
the inequahty (32) does not hold and there are no bound states of this type (type I). 
There is, however, another possibility for getting bound state. Since 

/(p + ire) = K + ix(p), /(p - in) = + ix(p), 

one gets only one real equation for real pi^2 if one puts Pi = Pi + m, P2 = P2 — ii^, 

X(pi) - X(P2) = 0. (33) 
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Noting that 

X(0) = X(7r) = 0, 

it is easy to see that there should be some value po at which x{p) has a maximum on the 
interval [0, tt] and the corresponding p'q — 2t:— po at which xip) has a minimum on the 
interval [tt, 27r]. As a matter of fact, po — There are no other extrema of x(p) on the 
interval (0, tt) . The presence of a maximum means that the equation 

xiP) = Xo 

has two distinct real roots if < Xo < X (^)) ^ < Pi < and < p2 < These 
roots serve also as nontrivial solution to the equation (33) and thus give the bound state of 
type II in which the wave function oscillates and decays exponentially as |ni — ^ oo. 
For Per < P < IT such a solution always exists. Similar solutions corresponding to 
— X (^-^) < Xo < can be found with any —tt < P < —Per- 

The above treatment is universal with respect to parameter k in the interval < 
K < 00. In the nearest-neighbor limit k ^ 00, the type II states with complex relative 
pseudomomentum and oscillating wave function disappear (Per- ^ ti" ) and the result 
coincides with the known one for the Bethe solution. 

3.2. Multimagnon scattering. After solving two-magnon problem, it is natural to 
try to find a way to describing scattering of M magnons with M > 3, i.e. find solution 
to the difference equation 

M 

J2 J2 V{np- s)tlj{ni,...np_i,s,np+i,...nM) 



M 

Y,V{np-n^) + J-^eM-Meo , (34) 

where n e Z^, the notation Z[„] is used for the variety Z — (ni, ...tim) and £0 = ^j^o ^U)- 
The exchange interaction V{j) is of hyperbolic form (5). 

The first attempt to solve (34) for M > 2 was made in [18] with the use of trial solution 
of the Bethe form and taking into account by semi-empirical way the corrections needed 
due to non-local form of exchange in (34). In this paper, the explicit solutions have been 
found for M=3,4 but the regular procedure of getting solution for higher values of M was 
not proved rigorously. The rigorous treatment of the solutions to (34) has been found 
later [21]. It is based on the analogy of the solution to (34) and corresponding solution to 
the quantum Calogero-Moser M-particle system with the same two-body potential and 
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specific value of tlie coupling constant in (3), determined by / = 1- This analogy is already 
seen in the form of two-magnon wave function (26) and holds for M = 3, 4 too. It was 
the motivation of the paper [21] to use this analogy in detail. 

The solution to M-particle system with hyperbolic potential and coupling constant 
with l—l is not simple too. The first integral representation for it has been obtained in 
[19] and more simple analytic form based on recurrence operator relation was given in 
[20]. I will follow [20, 21] in description of the M-magnon problem on an infinite lattice. 

Let us start from continuum model (3) with the interaction (5) and l—l. The solution 
can be written in the form 

) (x) = exp (^t tp,x)j ) (x) , (35) 
where ip^\x) is periodic in each xj, 

In [20], the explicit construction of the differential operator which intertwines (3) at (5) 
and / = 1 with the usual M- dimensional Laplasian has been proposed, and the functions 
of the type (35) have been represented in the form 

xf'\x)^DMe-K^\lY.P^^^A^ Dm ^ Q'm^-'Dm-i, (36) 



where 



ll...lm-l 

n 



d d 



2k COth K(Xj^ - Xn) 



+ Y: 2«:^ smh-'[^{xi^ - XiJ]Qil-''-'''^'-''-\ Qn = 1- (37) 



m—l 



=1 



This double recurrence scheme is very cumbersome because of presence of multiple dif- 
ferentiations but it allows one to reduce the construction of xi^^K^) to a much simple 
problem of solving the set of linear equations. Indeed, it follows from (36) and (37) that 
the function (p^{x) from (35) can be represented in the form 

<pI^\x) = i?({coth k{xj - Xk)}), (38) 

where i? is a polynomial in the variables {coth K,{xj — xi)}. As it can be seen from 
the structure of singularities in (3), the function (f^ {x) has a simple pole of the type 
[sinh K{xj — Xfc)]^^ at each hyperplane Xj — Xk = 0. As a consequence of (38), all the limits 
of if^^^x) as Xj ±oo, must be finite. Combining these properties with the periodicity 
of </?, one arrives at the following formula for the eigenfunctions of the Calogero- Moser 
operator: 

{M 1 M 

Y,[ip, - k{M - l)]x, n sinh-^ k{x, - x,)Sf\y), (39) 
/i=l I iJ.>y 
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where S^^\y) is a polynomial in = cxp(2Ka;^) in which the maximal power of each 
variable cannot exceed M — 1. Hence this polynomial can be represented in the form 



M 



(40) 



where is the hypercube in Z^, 

m e ^ < < M - 1, 

and dm{p) is the set of coefficients; it will be shown, however, that most of them 
vanish. The eigenvalue condition for the function (39) can be written in the form 



M r 

E 

13=1 L 



d 



d 



\^^dip ^ ''^'^^^ - M + IJ - m-'p0{M - 1) + (M - 1)(2M - l)/3 



M 

E 



yp-ypV dyp 

It can be satisfied if for each pair (/9, p) the polynomial 

d 



= 0. 



g{M) 



(41) 



yp 



dy0 dy^ 



yp^ + ^{pp - Pp) 



is divisible by {yp - yp). With the use of (40) this condition gives (M - 1)(2M - l)M^/2 
linear equations for the coefficients dm{p), 



y ] dmi...mp+n...mp—n...mM^P) 



mp -mp + 2n + —{pp - Pp) 



= 0. 



(42) 



The sum over n is finite due to restrictions to the indices of dm{p)- Substituting (40) gives 
also the set of equations 



M 



M r 



E IKM^-^pHE 



^ 2 2i / i 2M - 1 
2m^ H ppmp - 2m^ + -pp ; 



K 



Vm-1) 



M 

E 



yp + yp 



m/3 - rup + —(pp - Pp) 



} = 0. 



(43) 



p^p yp - yp 

After performing exphcit division by {yp — yp) in (43), one gets finally the second system 
of equations. The structure of the set dm{p) is specified by following propositions 
(for a sketch of proofs, see [21]). 

Proposition 3.1. S'^^^y) is a homogeneous polynomial of the degree M[M — l)/2. 

Proposition 3.2. The set of djn{p) can be chosen as depending on p and k only 
through combinations K^^{pp, —Pu)- 

Proposition 3.3. Let {P} be the following set of numbers {m^}: = Pfi — 1, 
where P is an arbitrary permutation of the permutation group ttm and I < p < M. The 
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nonvanishing dm{p) with coinciding values of {m^} arc expressed through d{pj{p). The 
latter are determined by the system (42) up to some normalization constant do, 



d{p}{p) =do]l 1 + ^(^^^-V - Pp-'u) 



(44) 



Proposition 3.4. Let (—1)^ be the parity of the permutation P. If xp(^+i) — — > 
oo, 1 < /X < M - 1, then 

hmxf H^)exp |^-^EP/3X^j = (-l)^2^^d{p-ij(p). (45) 

According to Proposition 3.3, the solutions to (42) must obey (43), and (43) has to 
be considered as a consequence of (42). Direct algebraic proof of this fact is still absent. 

The problem is now to solve the equations (42). It can be done exphcitly for M = 3, 4 
as follows: let [/Xi.../xm] be the permutation (1 — > /ii, ...M — > /im) and r^i, — i(2K)~^(p^ — 
p,^). Then, at M = 3 there are 6 coefficients of the d{py type which are caculated by the 
formula (44), 

douip) = c?o(l + ri2){l + ri3)(l + r23), c?io2(p) = do{l + r2i)(l + r23)(l + ^3), 

d2w{p) = do{l + r32)(l + r3i)(l + r2i), rfo2i(p) = do{l + ri3)(l + ri2)(l + r32), 
c?i2o(p) = do{l + r3i)(l + r32)(l + ri2), ci2oi(p) = 4(1 + ?^23)(1 + '^2i)(l + ^31). 
The only nonvanishing coefficient of another type is determined from (42): 

diu{p) = do{6- - rf^ - r^g). 

At M = 4, there are 24 coefficients of d{py type and other nonvanishing terms with coin- 
ciding values of indices can be arranged in three sets. The first two are given by elements 
with three coinciding indices and can be found from (42) by using known expressions for 
d{p} type, 

c^iii3(p) = do{l + ru){l + r24)(l + r34)(6 - - r^g - r^s), 

d222o{p) = do{l + r4i)(l + r42)(l + r43)(6 - rf^ - r^3 - rj^) 

and other elements of these sets dn^iip), di3u{p), dsniip) and ^2202 (p) , (^2022 (p) , (^0222 (p) 
can be obtained by the permutations [1243], [1342], [2341] of indices in these expressions. 
The remaining set consists of the coefficients with two pairs of coinciding indices, 

C^1122(p), d22u{p), C^2112(p), C?122l(p), d^uip), C?212l(p)- 

They may be determined by (42) with the use of known coefficients belonging to the first 
set, 

C?1113(p)(-2 + r34) + C?ii22(p)r34 + c?ii3i(p)(2 + r34) = 0. 
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and others come from the analogous equations arising after the permutations [3412], 
[3214], [4123], [1324], [4123] of the indices. 

These examples show that the solutions to (42) are crucial for determining the whole 
function Xp^^H^)- The question is now to see how these findings can be used for spin 
problem, i.e. the solution to the difference equation (34). The motivation is the striking 
similarity of the wave functions for M = 2. Guiding by it, I proposed the multimagnon 
wave functions similar to the functions like (35) with the structure (39), which are properly 
symmetrized combinations of them, 

ip{ni,...nM) ^ ]]_[sinhK{ni^-n^)]~^ (-l)^exp (i ^pp^^-A ) 



^ dmi...m„ 



ip) exp 



M 



A=l 



(46) 



where {d} is the set of unknown coefficients which might be determined from the M- 
magnon eigenequation if this Ansatz is correct. To verify the hypothesis (46), one has to 
calculate the left-hand side of the equation (34) with wave function of the form (46), 



M 



L{{n}) = K^Y^ [smh. ninp - s)] ^^/^(ni, n^j.i, s, n^+i, ...tim) 

/3=1 s€Z,„, 



M 



13=1 Pe-KM 



M 



W sinh K{n^ - n^) 



meDM 



where 



X exp <j ^ [ipp^ + K{2mp^ -M + l)]n^ \ W{ppp, rupfj, {n}), (47) 



W{p,m,{n}) = :j:x:2 



K 



M 



sinh ^ K{nx — s) 



,gZ[„, sinh k{s - np) ^^^^ 

xexp{[ip + K(2m- M+ l)]s}. (48) 

The sum (48) converges for all m G if p G C is restricted to |$5mp| < 2k. The explicit 
calculation of the sum (48) is based on the calculation of the function of a complex 
parameter x G C, 



exp(gs) 



M 



V- 

sinh^ k{s -np + x) 



[sinh n{nx — s — x)] ^, 



q = ip + n{2m- M + 1). 
As it follows from definition, this function is double quasiperiodic, 

Wg{x + inn-^) = exp[i7r(M - l)]Wg{x), Wg{x + 1) = exp{-q)Wg{x). 
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Hence it can be treated on the torus Ti = Z/Z + i7r/s:~^Z, and its only singularity on this 
torus is the double pole at a; = which arises from the terms with s = rii, ...um- The 
first three terms of its Laurent decomposition can be found directly from definition, 

Wg{x) = box'"^ + bix'^ + 62 + 0{x), 

M 

bo = exp(gn^) [sinh/«(nA - 

M 

bi — k{6o X] coth«;(n-y — n^) — ^ exp{qnp) 

M 

x[sinh«;(n/3 — Up) Y\_ sinh K{nx — np)]"^}, 

lM-11^ 
-- H ^ o 2^ coth(nT, - np) coth(n5 - np) 



K 



+ ^ sinh" 



E- 



exp(gnp 



^sinhfi;(n^-np) 



W [sinh n{nx - Up)] ^ 



X 



M 



+ W{p,m, {n}). 



coth — np) + ^ coth «;(n-y — Up) 

The next step consists in constructing the function Uq{x) with the same quasiperiodicity 
and singularity at x = by using the Weierstrass functions pi {x) , {x) and cti {x) defined 
on the torus Ti, 

Ug{x) = -A^ ^-exp{Sx) 

ai(x — r) 

x{p,{x) - p,{r) + A[(^{x + r) - ({x) - C(2r) + C(r)]}, 

where A, r, 6 and A are some constants and the term in braces is chosen as double periodic 
and having a zero at a; = r. Hence the only singularity of Uq{x) on Ti is double pole at 
X = for all values of r and A. 

Using the quasiperiodicity of the Weierstrass sigma function one gets 

cri(x + r + l) ,cri(x + r) ai(x + r + iiTK'^) .(Ti(x + r) 

— — ( = exp 277ir — , — — — ^ = exp 2?72r — r , 

ai[x — r + 1) ai{x — r) ai[x — r + nTK ai[x — r) 

where rji = 2Ci(l/2) and 772 = 2Ci(i7r/2K). Comparing these expressions with quasiperi- 
odicity of Wq{x), one finds two equations for r and S, 

2r)ir + 5 — —q, 2rj2r + mn'^S — m{M — 1). 

Their solution can be easily found with the use of the expression for q and Legendre 
relation rnK'^rji — r]2 — 27ri, 



m tp 
2"^4^ 



M 



1 + -rCi 

TT 



2^) 
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The Laurent decomposition of Uq{x) at a; = is obtained with the use of standard 
expansions of the Weierstrass functions, 

U,{x) = A[x-^ + (2Ci(r) + 5 - A)^-^ + ^(20 (r) + 5 - 2A)(2Ci(r) + 5) 

+A(2Ci(r)-Ci(2r))-pi(r)] + 0(x). 
The function Wq{x) — Uq{x) is analytic on Ti if A and A obey the conditions 

A^bo, A{2Ci{r) + S - A) = h. 

The only analytic function which is double quasiperiodic on the torus Ti is zero due to 
the Liouville theorem. Comparison of third terms in the decompositions of Wq{x) and 
Uq{x) gives the explicit expression of 62 in terms of bo,r,S and A, 

&2 = 6o[l/2(2Ci(r) + S - 2A)(2C. + S) + A(2Ci(r) - Ci(2r)) - pi(r)]. 

It allows one to find the explicit expression for the sum (48) in terms of p, m, {n}, 

{M 
— exp(qni3) [sinh K{nx — np)'^ 



where 



'M - 1) 1 ^ 

^- — h - ^ coth K{n^ — np) coth K,{n^ — np) 



(49) 



M 

+ [sinh K{n^ — ri^)]~^ — K~^f{r) ^ cothK(n-y — np) + K~'^e{r) 
^ exp(gnp) T-rr . , / xi-i 

coth K,{np — Up) + ^ coth K,{n^ — Up) — n^^ f{r) 

/» = Ci(2r)+5, 

^» = -y-^Pi(2r) + i/(r)^. 

It is worth noting that / and e are some polynomials in m. Indeed, it follows from the 
definition of r and b that 



m 



M -1 mC — 



where 



ip 
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By using quasiperiodicity of Ci{^) 

Ci(x + = Ci(^) + 2/C(l/2), 
one can represent the above functions as 

f{r) = f{p)-n{2m + l-M), 

where 

/W = C:(2r,) + *,= ^C,g)-C.(|). (50) 

Note that this function just coincides with the function (28) used for analysis of two- 
magnon scattering. The corresponding formula for i reads 



K 



e{r) = e{p) - n{2m + 1 - M)f{p) + — (2m + 1 - M)\ 



(51) 



where 



Now, according to (49-51) the left-hand side (47) of the eigenequation can be represented 
as follows, 

L(M) = Li(M) + L2(M) + L3(M), 



where 



M M ^^.2 

I]^(p/3) - — ixT ^ 



M 



M 

X ^ exp 



M 



[^Pp7 + K,{2mpy — M + l)]n^ 



X exp{[i(pp/3 ppp) 2hi{mpf3 mpp - M l)]np}[sinh «;(n/3 - Up)] ^ 



X 



M 



coth K{np — Up) + ^ coth «;(n-y — n^) — n ^fipp/s) + 2mpp — M + 1 



il sinh K(n.y — rig) 
X _[_[ 



M 



Lsiin}) = -K^ '[[[sinh K{n^-n^)] ^ ^(-1)^ J2 dmi...mMiP) 



M 



X exp < ^[ippj + K,{2mpj - M + 1)] 



.7=1 



(52) 
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M 

E 



^ K-'{2mpf, -M + l)f{ppp) + ^ ^ 

— y^[^~^/(pp/3) + M — 1 — 2mp/3] coth«;(n-y — n^j) 



+ ^ coth«;(n-y — n/3) coth«;(ny — n/3). (53) 

Now one can see that Li{{n}) exactly coincides with the right-hand side of the equation 
(34) if the M-magnon energy is chosen as 



M M 
/3=1 13=1 



2k J 



TT 



2aJ 



The problem consists in finding the conditions under which L2,3({n}) vanish. Consider 
at first the equation (52) and denote as Q the transposition /5 <-> p wich does not change 
all other indices from 1 to M. The sum over permutations in (52) can be written in the 
form 



M 



L2({n}) = -K^ Yl [sinh«;(n^ - n^)] ^ ^ c?^i...mM(p) 



M 



where 



F{P) = exp 



M 



J2 {Wpi + «(2mp^ - M + 1)) 



imp 



X exp{[i(pp^ + ppp + 2K{mpi3 + rupp - M + l)]^^} 

_i it sinli/v(//-, - 113) 

X smh «;(n;3 - Up) M . , , ^ 

^^^pS\nhK{n^-np) 

x-[2mp^ — K~^f{pi3) + cothft;(n^ — n^) + ^ coth«;(n^ — n^) — M + 1]. 

Note that the only difference of F{PQ) and F{P) is in first two terms in last brackets. 
This allows one to rewrite the last formula as 

L2{{n})^-K\X{[smhK{n^-n,)]-' ^ (-1)^ 

P'>f PSttm 



X ^ exp 

P¥^P 



^ [ipp^ + K(2mp^ - M + 1)] 



xsinh U{np-np) [] sll^wl"' _ 

2{M-1) 

X E E exp{[i(pp^+ppp) + 2K(s-M + l)]np} 

{m^}eDM ,'r^PI3,Pp s=0 
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X [il^ \s i\/f +1|] ^ ] ^ ] ...mpfi+n...mpp— n...mM 



X 



The comparison of the last sum with (42) shows that it vanishes if 



(54) 



where d^m}{ifip)) is an arbitrary solution to the system (20) with replaced by i/(p^), 
1<H<M. 

The only problem is now transformation of Ladn}). Taking into account the formula 



M 



^ coth K(nj - np) coth K{n^ - np) = -M{M - 1) (M - 2) 

and symmetrizing over /3,7 in (53), one finds 

-^3({^}) = -1^^ n [sinh«;(n^ - n^) 



X H (-l)^exp 



M 



^[ip^ - k{M - l)]np-iy 

7=1 



R{P,{n}), 



where 



M 



R{P: {^}) = Yl drm...mM 

meDM 



X 



M 

E 



l(M-l-2m,)^ + ^^ 



(p) exp \ 2k ^ np-i^m^j 

-K-^f{pp){2mp-M + l) 



- Y.\P^l3 - '^1 - - /(Pt))] cothK(np-i;3 - np-i^) 

Upon introducing the notation exp(2«;np-i^) = at fixed P, one finds 



M 



m. 

*'mi...mMyP) I 11 y7 

mGZ)Af \7=1 



7 



M r 



2M-1 



3 ; 



) {M - 1) 



M 



-E 



VlB + y-y 



[mp -ruj- {2k) \f{pp) - f{Pj))] 



Now it is quite easy to see that replacing dm^...mMiP) ^ dm^.-mMiP)^ ^f{P^l) ^ Pa* 
the right-hand side just gives the left-hand side of equation (43) and must vanish for 
all y e if the set d{rn} solves the equation (42), i.e. the function x^^^ satisfies the 
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Calogcro-Moser eigenequation. Hence both L2,3({n}) vanish under the conditions (54) 
and the Ansatz (46) satisfies the eigenvalue problem (34). 

These lengthy calculations lead to the simple receipt: to get a solution to (34), one 
needs to change the p dependence of the perodic part of the solution to hyperbolic 
Calogero-Moser quantum problem as {p — > ifip)}- The asymptotic behaviour of the 
M— magnon wave function ip{ni, ...um) (46)as «; — > oo or [n^u — nu\ — > oo can be found 
with the use of Proposition 3.4. In the former case one obtains the usual Bethe Ansatz 
[4,5] as a conseguence of (45) and the relation 

The generalized Bethe Ansatz appears at finite n when the distances between the positions 
of down spins tend to infinity as np(A+i) — npx — > +oo, 1 < X < M — 1, 

/ M 

t/'(ni, ...nju) = V^o (-l)^^exp U ^pqa'T'A 

QSttm \ A=l 

xn 1-^[/(PQPm)-/(PQp.)] , (55) 

where f{p) is given by the formula (50). The asymptotic form (55) will be used in the 
next section within the asymptotic Bethe Ansatz scheme of calculations of the properties 
of the antiferromagnetic ground state of the model. 

According to (55), the multimagnon scattering matrix is factorized as it should be for 
integrable models. There is a possibility for existence of multimagnon bound complexes 
for which some terms in asymptotic expansion (52) vanish. Such a situation does not take 
place for usual quantum Calogero-Moser systems with hyperbolic interaction where the 
two-body potential is repulsive. 

4. Periodic boundary conditions and Bethe-Ansatz equations 

Imposing periodic boundary conditions (with period A^) for the spin chains with inverse 
square hyperbolic interaction leads to the elliptic form of exchange (12). These conditions 
allows one to treat correctly also the important case of antiferromagnetic case which 
corresponds to the positive sign of coupling constant J in (12). 

The spectrum of one-magnon excitations over the ferromagnetic ground state is now 
discrete and can be calculated via Fourier transform of the eUiptic exchange [17]. Through- 
out this section, the notation u; — in /k, will be used for the second period of the Weierstrass 
functions. As in the previous Section, I will consider at first the case M—2 which allows 
more detailed description. 
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4.1. Two-magnon scattering. As in the case of infinite lattice, the problem consists 
in finding two-magnon wave function defined by 



1^ ^ V(^l,^2)s„iSnJO >, 

where j'^ > is an eigenvector of the Hamiltonian and |0 > is the "vacuum" vector with all 
spins up. The corresponding two-particle problem is now the Lame equation, and well- 
known Hermite form of its solution allows one to guess the Ansatz for the wave function 
in the form 

exp[i(pini + P2n2)]crN{ni -n2 + -f)+ e^[i{pin2 + P2ni)]aN{ni + 713-7) 

ip[ni,n2) — r . 

(7Ar(ni - 712) 

Since ip should be periodic in each argument, the parameters pi,2 are expressed through 
the phase 7, 

PiN - irji'j = 27rZi, P2N + irji'j = 277/2, (56) 

where 771 = 2C,n{N /2),r]2 = 2Cjv(a;/2) and li,l2 e Z. 

The solution to the eigenequation is now based on the formula 

aN{k-l + -f + z) (Tat (/- 7)0-1 + To ^ 



J2 pN{k + z) — — . exp(iaA:) 



«7 ) 



X exp 1^ [C7v(iV/2)Ci(o;/2)7 + tCN{oo/2)a] 

Cn{1 - 7) - Civ(0 exp(m/)ajv(7)ajv(0 
+ 2 2M/-7) ^'^^^ 

+^(Civ(A^/2)Ci(a;/2)7 + iCiv(a;/2)a)]}, 
where I e Z and a and 7 are connected by 

expfiaAT + 27Cjv(A^/2)] = 1, r^,^ = -(47r)-i[Q;u; - i7Civ(c^/2)]. 

The two-magnon energy is given by 

£2(Pl,P2, 7) = •>^{l/4[/(Pl, 7) + /(P2, -7)]^ + £o(Pl, 7) + £o(P2, -7) + P{l)}, 

where 

/b' 7) = Ci (^^^^^^) + (^^)"' [^2Ci(l/2) + ^pCl(^/2)] (57) 
and pi,2 and 7 are constrained by 

/(pi,7)-/(P2,-7)-2Civ(7) = 0. (57) 
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With the use of (56), (57) and direct computation it is possible to show that 

N 

S+ Yl i'{ni,n2)s-^s-JO>^0, 

i.e. these states have the total spin S — — N/2 — 2. 

It is natural to ask of how many solutions do the equations (56), (57) have. The 
completeness of the set of these solutions means that their number should be equal N{N — 
3)/2 since in two-magnon sector there are N solutions with ■0(^1; ^2) = '0i(^i) + '0i(^2)- 
Is it possible to evaluate the number of solutions to (56), (57) analytically? The answer is 
positive [22] . The sketch of the proof is as follows. The constraint (57) can be rewritten 
as 

+^[Civ(c^/2) - iV-^Ci(^/2)] - 2Civ(7) = 0. 

At fixed Zi,2, it is a transcendental equation for 7. 

Let now A be the manifold which consists of various sets {Zi_2 G Z, 7 e C} and call 
two sets {Zi, ^2, 7}, {I'll 7'} £ ^ equivalent if the corresponding wave functions coincide 
up to normalization factor. With the use of (56) and quasiperiodicity of sigma functions, 
one finds that the manifold A is equivalent to its submanifold Aq defined by the relations 

< /l < iV - 1, /2 = 0, 7 e T^N,Nu.- 

Let {A} be a variety of nonequivalcnt sets within Aq. The question now is: how many 
nonequivalent sets obeying 0(7) = are in Aq? To answer it, let us note that the 
function 0(7) is double periodic with periods N and Nuj and there is the relation 
between C, functions of periods {N,u!) and {N,Nu!), 

2x 

C^{x) = Cix) + J2 [C(^ + M - CiM] + -[Cn{uj/2) - C(iVa;/2)], 
3=1 ^ 

where ({x) is the zeta function defined on the torus T^^nui- With the use of scaling 
relation (i{N~^{x) = N((x) one can rewrite the constraint ^^^^0(7) = in the form 

N-i /Nu\ 
-2 E C(7 - M + mCil - hcu) + C(7)] + 2C — (/i - + 1) = 0. (58) 

j=0 V ^ / 

It is easy to see that at A?^ > 2 there are A^ simple poles of the left-hand side of equation 
(58) located at 7 = ju and this function is elliptic. Then there should be just N roots 
of equation (58) within TN,Nio ^'^d, at first sight, {A} consists of A^^ elements. However, 
some sets with different roots of (58) may be equivalent. In fact, one can see that if 
7o £ Tjv.ATw is a root, then 

7o = -70 + + A^sign(Ke7o) + - sign(/i|cc;| - '^m'jo)] 
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is also a root of (58). Moreover, all the solutions to the equation 7q = 70 are the roots of 
(58). The sets of these solutions are different for N, li even or odd. There are four cases. If 
both N and u; are even, there are only two these roots, (iV + /ia;)/2 and (A^ + /ia; + A^u;)/2. 
If both N and h are odd, the additional root liUj/2 is present. As N is odd and h even, 
the additional root is {N + li)u/2. And in the case of even N and odd I, one has four 
such a roots since both liUj/2 and {N + li)uj/2 obey the equation (58). 

All these exphcit roots are combinations of half-integer periods of the torus Tjv.atw. 
In this case the wave function can be simplified and it turns out that ip{ni,n2) vanishes 
identically for all exphcit roots listed above. 

The number of all nontrivial and nonequivalent sets {^1,0,7} in the variety {A} can 
be now easily counted. At even N, there are A^/2 even {li} with {N/2) — 1 nonequivalent 
roots and A^/2 odd {h} with (A^/2) - 2 ones. At odd N, there are (N + l)/2 even {h} 
and {N — l)/2 odd {li} with {N — 3)/2 nonequivalent roots in both cases. Hence the 
total number of elements in the variety {A} equals N{N — 3)/2 as it should be, and the 
nonequivalent solutions to (58) provide complete description of nontrivial two-magnon 
states. It would be of interest to investigate, in the limit of large N, the distribution of 
nontrivial roots within the torus Tjv tv^. 

4.2. Multimagnon states. As in preceding Section, one has to investigate first 
the solutions to usual quantum Calogero-Moser problem with coupling constant / = 1 in 
(3) and elliptic two-body potential. This problem at M > 2 was attacked first in [24] 
where the general statements on the structure of many-particle wave function have been 
proved and explicit result for M = 3 has been obtained. The problem of arbitrary / 
and M = 3 has been considered in [25] and soon the analytic expression for arbitrary M 
has been also found [39] in the process of solving the elliptic Knizhnik-Zamolodchikov- 
Bernard equation. Unfortunately, its form turned out to be so complicated that no explicit 
calculation were possible for multimagnon wave functions. At M = 3, the 3-magnon wave 
function has been found explicitly in [31] but the calculations were very lengthy and it has 
been not seen how to generalize the method for M > 3. The way to the solution of the 
M-magnon problem which does not refer to explicit form of the solution to M-particle 
problem has been found later [30,40]. Before describing it, it will be of use to formulate 
basic facts about the wave functions of the continuous M-particle problem for elliptic 
two-body interaction [24]. 

Since p{x) is double periodic, it is easy to see that the corresponding M- particle 
Hamiltonian (3) commutes with 2M shift operators Qaj — explujad/dxj), where 001^2 are 
two periods of p{x). Let x^K^i^ ■■■Xm) be their common eigenvector, 

2 2 M 2 

X^^\x, + E t^oo^, -XM + E ^M^o) = exp(zE Epi''^4"^)x^'H^i, -^m), 

a=l a=l j=l a=l 

where /j"'' G Z. Hence x''^\x) can be treated on the M-dimensional torus Tm = 
C/Zlui + Zu!2 with quasiperiodic boundary conditions. The structure of singularities 
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of the Hamiltonian (3) in this torus shows that 

^(p) is 

analytic except of all hypersurfaces 
Ljk defined by the equalities Xj = x^, 1 < j < k < M. On each Lj^, x^'^^ has a simple 
pole. Let be a class of functions with these properties. 

Proposition 4.1. The class is a functional manifold of dimension 2M — 1 + 
rpj^g parameters {p^^} are not independent but restricted by the linear rela- 
tion J2jLi{Pi''<^2 — P2^i^i) — Zo^i + Zu2. The manifold can be described as an union 
of the (2M — l)-parametric family of hnear spaces with dimensions M^~^ with the basic 
vectors parametrized by {p^-*}. 

Proposition 4.2. The co-ordinate system on can be chosen by such a way that 
all its elements are expressed through the Riemann theta functions of genus 1 or usual 
Weierstrass sigma functions. 

The sketch of the proofs can be found in [24]. The explicit expressions for x^'' can be 
also found in [24] for M = 3 and in [39] for arbitrary M. The amazing fact is that the 
treatment of M-magnon problem can be done without use of these explicit expressions. 
Let us choose the exchange in the form 

2 



h{j) = J 





7r\2 r 


(- 


sm — 


Vtt 


CO J . 



Pn{j) + -Cn (I) 



so as to reproduce correctly the inverse square hyperbolic form of the Section 3 in the 
thermodynamic limit N oo. The second period of the Weierstrass function p^v is 
uj = i-n / K. The eigenproblem is decomposed into the problems with M down spins due 
to rotation invariance and the eigenvectors > are given by 

N M 



|^(^)>= Yl ^M(ni..nM) n ^nJO >, 



ni..nM 



f3=l 



where |0 >= | ft ... t> is the ferromagnetic ground state with all spins up and the 
summation is taken over all combinations of integers {n} < N such that Y\^^j,{nn — 
rii,) ^ 0. The function ipM is completely symmetric in its arguments and obeys lattice 
Schrodinger equation 



JV M 
s^ni,..nM /3=1 



s)'0M(m, -n^-i, s, np+i, ..um) 



+ 



M 



and the eigenvalues of the Hamiltonian are given by 



(59) 



'UJ . ttV 
J ( — sm — 



H 



2M(2M -l)-N 



UJ 



Let Xm be the special solution to the continuum quantum many-particle problem 



1 



M 



-2^ 



M 



. f, 2 +Y1 Pn{xi3 - Xx) - Em(p) 
/3=1 /3^A 



Xm(xi,..xm) =0, 
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which is specified up to some normaUzation factor by particle pseudomomenta (pi, ..-Pm) 
and obeys the quasiperiodicity conditions 

Xm (xi, ..xp + N, ..xm) = exp(ip^Ar)xM (xi, ..xm), 1 < /5 < M, 

Xm ..X/s + to, ..Xm) = exp(27rig^(p) + ippuj)xM i^i, ..Xm), < ^e{qp) < 1. (60) 

As will be seen later, the set {qisip)} is completely determined by {p}. 

The connection of Xm with multimagnon wave function is given by the Ansatz 



M 



/ M \ 

exp \ -i Y^p^nA Xnini, ..um), 



(61) 



where 



p^^p^- 2t:N e Z. 

The last condition is just the condition of periodicity of ipM- The problem now consists in 
calculation of the left-hand side of the lattice Schrodinger equation (59), but before doing 
this let us mention that Xm(p) has the singularities in the form of simple poles and can 
be presented in the form 



M 



(p) F^p\xi,..xm) ^, \ T[ ( ^ 

U[Xi,..Xm) a<fi 



(62) 



where (Jn{x) is the Weierstrass sigma function on the torus T/y. By definition, the only 
simple zero of (JNix) on T^r is located at a; = 0. Thus [G{xi^ ..Xm)]~'^ absorbs all the 
singularities of Xm ^^^^ ^^e hypersurfaces x^ = Xp. The numerator F*^^) in (62) should be 
analytic on (Tjv)^^. It obeys the equation 



M 

a=l "-^a 



+ 



M 



M 



2Em{p) - ^ X) {pN{Xa - Xf}) - Clf{Xa - Xp)) 



pip) 



C^iXa - Xp) 

V 



dp(p) dF^P^ 



dXr 



]. 



The left-hand side of this equation is regular as x^ — > x,^. Hence F^^ must obey the 
condition 



d 



d 



^dxn dxt,J 



]f^p\xi,.. 



Xm)\x^=x^ 







(63) 



for any pair (//, i^). Let us now show that the properties (60), (62,63) allow to validate 
the ansatz (61) for ■^m- Substitution of (61) to (59) yields 



M 



PeTTM 1/3=1 



M 



PN{np/3 - npj) - Em 



Vui'Tipu -npM) 



0, 
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where 

N 

Sf3{npi, ..npu) ^ X] pN{npf) - s)Q^i^^(p^{npi,..npM) (64) 

ST^npi,..npM 

" (s) 

and the operator Qjj replaces /3th argument of the function of M variables to s. 

The calculation of the sum (64) is based on introducing the function of complex 
variable x 

N 

Wi,^\x) = X pN{npp -s- x)Q'j^^'''^ip^^\npi, ..upm)- 

s=l 

As a consequence of (60) it obeys the relations 

Wi,^\x + 1) = W^/\x), Wi.^\x + 00)^ ex.p{27riq^{p))wi.^\x), (65) 

where ^ 

qp{p) = 5/3 (p) + 

The only singularity of Wp^^ on the torus Ti = C/Z + Zu is located at the point x — 0. 
It arises from the terms in (64) with s — npi, ..npM- Hence the Laurent decomposition 
of Wp'' near x — Q has the form 

W'j,^\x) = w_2X~^ + w_ix-^ +wo + 0{x). (66) 

With the use of (62), one can find the explicit expressions for w^i in the form 

W-2 = </?M (npi, ..Upm) 

w-i = ^ <y^M (^Pi> -npu) 

onpp 

+ (-l)^G'"^(ni, .Mm) XI Tj3x{npi, ..npM)Qp'''^^ exp -i^pyUpy j F^\npi, ..upm) 

1 

Wo = Sp(npi, ..Upm) + ^'dr^^^^^'^^^' ■■'^pm) + (-l)^G'"^(ni, ..um) 
X X Tfjx{npi, --Upm) Uf^xinpi, ■■npM)Qp'''^^ + pNiupp - npxjdQ^^''^^ 

X^0 

X exp I —i X PylT-Py 

Tpxinpi, ..Upm) = crN{npx - npp) J| ^^(J^Pp — ^}Pli^ 

p+(i,x ^N\npp - npx) 

Ui3x{npi, ..Upm) = p'n{^px - ^ph) - pN{np0 - npx) X (.N{npp - npx), 



where 
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(— l)'^ is the parity of the permutation P and the action of the operator dQ^p^^^^ on the 
function F of M variables is defined as 

d 

Note now that the expression for the function w'p\x) obeying the relations (65) and (66) 
can be written analytically without any further freedom, 

(7"l(r/3 - x) 

x[Ci(^ - rp) - Q,{x) + Ci(r^) - Ci(2r^)]}. 
The parameters a^,r^ are chosen as to satisfy the conditions (65), 

0/3 = 2g/3(p)Ci(l/2), r0 = ~qp{p)- 

By expanding the above form of W^^ in powers of x one can find Wq in terms of 
W-2i qp and obtain the explicit expression for (S^(npi, ..upm) ■ After long but straight- 
forward calculations the equation (59) can be recast in the form 



E 



ip^^\npi, ..Upm) 



= \g ^{ni, .Mm) Y. E [^/3a(?T'Pi, -^ipm) + Zxf^iripi, ..upm)] , (67) 

where 

//3(p) = 2g>(p)Ci(l/2)-Ci(g/3(p)), 

and Zpx{npi, ..npu) is defined by the relation 

Zi3\{npi, ..Upm) = Tpxinpi, ..upm) Upxinpi, ..npM)Qp'''^^ + pN{npx - npp) 

x{dQ^^''^ - U{p)Qt"% exp [-lYp.npA F^\np,, ..upm). 

One observes with the use of the definition (61) of that each term of the left-hand 
side of (67) has the same structure as the left-hand side of the many-particle Schrodinger 
equation and vanishes if £m and f^ip) are chosen as 

ffsip) ^ -ipfs, /3 = 1,..M, (68a) 

M 

£M^^M{p) + J2^piP)- i^^b) 
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It remains to prove that that the right-hand side of (67) also vanishes. It can be done by 
using the observation that the sum over permutations in it can be simply recast in the 
form 

J2 (-1)^ ■■npM) - ZxfsiripRi, ..uprm)], 

where it! is the transposition (/3 <-> A) which leaves other numbers from 1 to M unchanged. 
Taking into account the fefinition of Z, one finds 

Zf3\{npi, ..npu) - ZxpinpRi, ..uprm) = Tfjx{npi, ..npM)pN{npx - npp) 



X exp 



M 

{Pl3+Px)npx+ Pp'^Pp 



Now it is clearly seen that the last factor vanishes due to the condition (63). 

The relations (68a-b) for the spectrum are still not complete since the dependence of 
{q} on {p} is not known on this stage. This completion can be done only by further anal- 
ysis of the properties of solving M-particle Schrodinger equation. In [39] the explicit 
form of Xm{x) has been found in the process of solving the Knizhnik-Zamolodchikov- 
Bernard equations. In suitable notations, it reads 



M 

/3=1 seiTr, ^ 



where m — M{M — l)/2, c is non-decreasing function c : {l,..,m} — > {1,..,M — 1} 
such that |c~^{j}| — M — j, l{s) is an integer which is defined for the permutation s by 
the relation Xc{s{i))+id/dxc{s{i))--Xc{s{m))+id/dxc{s{m))xf = 1{s){xi...xm), {t} is a set of m 
complex parameters obeying m relations [39] 

E P(h -ti)-2 J2 Pih - ^i) + M^cjMh) = ^(PcO) - Pc(j)+i), (70) 

l:\cil)-c(j)\=l lU^j,c{l)=cij) 

p{t) = CAr(t) - ^CNiN/2)t, 

and 

a^{t) = exp{{2/N)CN{N/2)wt) ^^^"^ " ^' 



The main advantage of the explicit form of x function is that it allows to find the second 
set of relations between the Bloch factors {p}, {q}. It is easy to see that {p}'s in the 
definitions (60) and (69) are the same. The problem consists in calculation of {q}. To do 
this, it is not necessary to analyze each term in the sum over permutations in (69) since 
all of them must have the same Bloch factors. It is convinient to choose the term which 
corresponds to the permutation 

So ■ So{j) = m + 1 - j, J = 1, ..m. 
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After some algebra, one finds that this permutation gives nontrivial contribution to (69) 
with 1{sq) = M\[M — 1)!...2!. Moreover, with the use of explicit form of the color function 
one finds 

c(so(l)) = M-q if q{q-l)/2 + l <l<q{q +l)/2. 

Now the problem of calculation of the second Bloch factors reduces to some long and 
tedious, but in fact simple calculations of the product of factors which various a functions 
acquire under changing arguments of x function to the quasiperiod a;. The final result is 
surprisingly simple, 

qpip)^N-'( ti- E 1</3<M-1, (71) 

\/:c(0=/3 /:c(0=/3-l / 

with the first and second term being omitted ior P — M and P — 1. 

The equations (71), together with (68a) and (70), form a closed set for finding Bloch 
factors {p}, {q} at given integers {Ip} e Z/MZ and determining the eigenvalues of the 
spin Hamiltonian completely. The corresponding eigenvalue of the continuum M-particle 
operator is given by [39] 

13= 



N 



-Civ(Y) + Ep?/2 
^ ^ ^ 0=1 



E(2'^c(fe),c(0-^(^ifc - h) — ^\c{k)-c{}.)\,iF{tk — k)) - E ^i^k) 

k<l c(fe)=l 



where 



F{t) = -p^it) + {(Nit) - 2/NC^{N/2)Y + 4/NCn{N/2). 

It is worth noting that for real calculation of the eigenvalues one has to solve the Bethe- 
type equations (68a), (70), (71) at first. It is not clear how to treat properly this huge 
system of highly transcendental equations even in the limit iV — > cxd. In this limit, there 
is a procedure known as asymptotic Bethe Ansatz (ABA) which consists in imposing 
periodic boundary condition on the asymptotics of the wave functions for infinite lattice 
[41]. It will be used in the next subsection for obtaining some results on antiferromagnetic 
ground state. 

4.3. ABA results for large N. In this subsection, the ABA hypothesis (still un- 
proved) will be used for description of some properties of the spin chain with the exchange 

sinh^ K 



Kj) 



(72) 



sinh^ K{j — k) ' 

which corresponds to J = — (sinhK/K,)^ in (23) (the antiferromagnetic regime) at large 
but finite A^. Note that in the nearest-neighbor limit k ^ oo one can decompose (23) 
with the exchange (72) as 



E(^i^.+i - 1) + \^~"^ E(^i^.>2 - 1) + o(e- 

3 3 



2k 



)■ 
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Hence one can write the ground-state energy per site as 

e = ^ < a,,T,+i - 1 > +^e-'^ < - 1 > +o(e-2»), (73) 

where < > means average on the vacuum state of the Hamiltonian. Fortunately, in the 
first order approximation one can replace this state to the vacuum state of nonperturbcd 
Hamiltonian with the interaction of nearest-neighbor spins, Hq — \ j2j{^j^j+i ~ !)• It 
gives an opportunity to find the second- neighbor correlator by calculating (73) 

explicitly. 

The scheme of ABA is based on asymptotic expression of the wave function with M 
down spoins for infinite chain in the region rii <^ n2... <C um, which has been described 
in Section 3, 

/ M \ (i ^ \ 

V^(ni, ..Um) oc exp U ^ kpaUa exp - ^ x{PPa,Ppp) , 

PGttm V a=l / \ a<l3 J 

where the first sum is taken over all permutations from the group hm, {Pa} is the set of 
pseudomomenta and x{Pa,Pi3) is the two-magnon phase shift defined by the relations 

.X(Pa,Pp) . ^ . ^ 

cot = ip{pa) - ^{Pp): 

To consider the chains of finite length N in the thermodynamic limit N ^ oo,we adopt the 
main hypothesis of ABA, i.e. imposing periodic boundary conditions on the asymptotic 
form of the wave function. Taking ip{n2, ..nM,ni + N) = ip{ni, --nM) and calculating the 
both sides with the use of the above formula for asymptotics results in the ABA equations 

/ M \ 

ex^{ipaN) = exp U ^ xiPa^Pa) , a = 1, .M. (74) 

V Pl^oc J 

The energy of corresponding configuration is given by 

-Em = £ ^^"^ " {cos{kan) - 1). 
a^in^smh an 

For investigating the antiferromagnetic vacuum of the model one should take N even, 
M = N/2. Taking logarithms of both sides of (74) and choosing the proper branches, one 
arrives at ^ 

1^ ^ ~ ^ J arctan[(^(pe«) - 'fiipp)], 

where {Q} is the set of (half) integers. For antiferromagnetic ground state, one assumes as 
usually that these numbers form uniform string from —Qmax to Qmaxi Qmax = 

iV/4-1/2 
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without holes. After introducing rapidity variable A by the relation A = '^{k) and the 
function //(A) via the relation n — k = /i(A), the ABA equations (74) can be written as 
[38] 

QjN^ZiK), (74a) 

where 

Z(A) = (27r)-V(A) - 4? E arctan(A - A^). 

Following [38], let us go to continuous variable x — Qa/N in the hmit — > oo and 
introduce the root density crjv(A) by the relation (Tjv(A) = dx/dX. Differentiating both 
sides of (74a) with respect to A, one arrives at the following equation in the limit A" — > oo 

/CO 
A{X-X')a^{X')dX', (746) 
-oo 

where A{X) — [k{1 + A^)]~^. The energy per site can be written from as 

Coo = lim N-'En/2 = r e{p{X))ao.{X)dX, (75) 

where 

/ /x\\ r, • 1 2 cosnj9(A) - 1 
e{p{X)) = 2 smh" n J2 ■ 

smh KTi 

The solution to (74b) can be found via Fourier transform, 

/oo p^^^rlh poo 

Substituting it into (75) yields 

dXe{p{X)) / dk-——^^ / i,'{T)e-^'-dT. 



Upon choosing variables as A = (fi(p), ii'{T)dT — —dp' and changing the order of integration 
(it is allowed since the integral over r vanishes sufficiently fast as |A;| — > oo), one arrives 
at the following formula for an energy per site, 

/oo fJU /'27r f2n 

y—^J dpe{p)^\p)e^'^^^^ dp'e'^'^^^K (75a) 
-oo 1 + e 1*^1 Jo Jo 

where the functions e(p) and (p{p) are determined as above. Unfortunately, the integrals 
in (75b) cannot be evaluated analytically; however, one can find as k — > oo that 

(^(p) = Icot^ + 2e"^''sinp + 0(6"^"), 

e{p) = 2(cosp - 1) + 2e-2''(cos2p - 1) + oie'^"). 

Upon substituting these expressions into (75b), the inner integrals are calculated analyt- 
ically up to the order e~^'^ and final result for second-neighbour correlator in the model 
with nearest-neighbor exchange reads 

< ajaj+2 >= 1 - 161n2 + 9C(3), 
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where ( is the Riemann zeta function which appears in the right-hand side due to the 
formula /q°^ = 3/2((3). This result coincides exactly with the expression given by 
Takahashi [37] who considered the hmit of infinite one-site repulsion in the half-filled 
Hubbard model. 

Another ABA result is the calculation of central charge c of underlying conformal field 
model [23]. It is given by the formula for finite- A?^ correction to the energy of antiferro- 
magnetic ground state 

TTC 

where ^ is the velocity of the lowest-lying elementary excitations. The value of Ae^ 
can be calculated via the equations (74a) where the values of the order N^'^ should be 
carefully taken into account. I would like to mention only the final result of rather long 
calculations [23], 

AeAr = -(12iV2)-Voc + 0(Ar-3), 

0OO = 27r« lim V; ^^ ^ ■ 

The energy and momentum of elementary excitations over antiferromagnetic vacuum can 
be also calculated on the base of (74a) under an assumption that this excitation corre- 
sponds to presence of a hole in the sequence of numbers {Q}- These calculations result 
in the formula ^ = (27r)~Voo which gives the value of the central charge c = 1 as in the 
case of usual nearest-neighbor chain. 



5. Inhomogeneous lattices. 



It is generally believed that more general dynamical Calogero-Moser systems describing 
particles with internal degrees of freedom are integrable. The motion of particles can be 
eliminated by arranging them into classical equilibrium positions. By this way, the first 
model of inhomogeneous chain [27] has been obtained where spin interaction was given by 
inverse squares of distance between them and spins were located on equilibrium positions 
of particles with rational two-body interaction in the field with a harmonic potential. As 
for inverse hyperbolic square exchange, the integrability of the corresponding models is 
still questionable. Anyway, there are many indications to this fact as it will be shown 
later. 

The integrability of classical Calogero-Moser systems in some external fields has been 
considered in [35]. It was shown there that the Hamiltonians (3) with interaction (5) 
(with K = 1 as it can be removed by scaling transformation) are still integrable if the 
external field with the potential 

W{x) = cosh(4a;) + 2(3 cosh(2a;) + 27 sinh(2a;) (76) 
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is added. As for spin chains, the Hamiltonian is still given by 

N 

H = J2hjkPjk, (77) 

j<k 

where {Pjk} is any representation of the symmetric group ttn, hjk = sinh~^(a;j — Xk) and 
{xj} are the coordinates of classical particles at equilibrium obeying the equations 

-2j2hjkCjk + W'ixj)^0, (78) 

where 

Cjk — coth(xj — Xk). 

The first question is to construct the Lax pair for these systems. Consider the following 
Ansatz of {2N x 2N) matrices (L, M) with entries 

^11 ^ _jJ22 ^ ^12 ^ ^21 ^ ^ ^ ^11 ^ ^22 ^Mo + m, M^^ = M^^ = 0, 

where Lq and Mq is the standard Lax pair for the systems without external field, 

N 

(-f'o)jfc = (1 ~ Sjk)CjkPjk, (M))ife = (1 ~ Sjk)hjkPjk — Sjk hjgPjs 

and 0, p and m are (A^ x A^) matrices with entries 

(V')jfc = Cizj)Sjk, (pjk = (pizj)Sjk, {m)jk = IJ'izj)Sjk, {p)jk = (1 - Sjk)Pjk, 

where 2;^ = exp(2a;j). The Lax relation [H, L] — [L, M] is equivalent to the set of functional 
equations 

Cjk[l^{Zj) - IJ,{Zk)] + [ip{Zj) + ip{Zk)] = 0, 

Cjk[<p{zj) + (p{zk)] + hjk[^{zj) - ^{zk)] + ii{zj) - ii{zk) = 0. 
The general solution to this set is given in [36] , 

ll{z) ^ lIlZ + ll2Z''^ , (p{z) ^ -/IiZ + II2Z~^, C{z) ^ IIiZ + II2Z~^ 

The potential of an external field reads 

W{z) = 2[filz^ + filz-^ + (27 - l/2){fi^{z) + i^2Z-% 

It contains three free parameters as (76). For the special case of the external Morse 
potential {fi2 = 0) the matrix M obeys also the condition Y.'j=i ^jk = 0, which guarantees 
that the integrals of motion can be constructed as {}2f^{L'^)jk}- In other cases, the 
existence of the Lax pair does not imply integrability immediately. 
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The extra integrals of motion should be some polynomials in the permutations as it 
takes place for usual lattice spin models [32]. It turns out that minimal degree of this 
polynomial is now equals 3 and the operator 

commutes with if if F is a solution of functional equation 

g{xj, Xk) + g{xk, xi) + g{xu Xj) = 0, 

where 

gixj,Xk) = 2hjkiF{xj) - F{xk)) + CjkiW'{xj) + W'{xk)). 

The solution is given by the relation g(^Xj,Xk) = G{xj) — G{xk) and functional equation 
for the potential 

c^k{W'{x^) + W'{xk)) - 2h^k{W{x^) - W{xk)) = G{x^) - G{xk). 

Its general solution just gives the form (76) which supports the hypothesis of complete 
intcgrability of this class of models. 

To construct the explicit eigenvalues of the corresponding spin Hamiltonians, one needs 
more knowledge about the solutions to equilibrium equations (78). It can be easily done 
for special case of the Morse potential W{x) = 2r^(exp(4x) — 2exp(2a;)), where these 
equations have the form [26] 

-E¥^^ + -'(-.-1) = 0, (78a) 

where the variable z — exp(2x) is introduced. Following the observation of Calogero 
[28], one can assume that the roots {zj} of (78a) are given by roots of some polynomial 
Pn{z) — Y[jLi{z — Zj) obeying the second-order differential equation. In the case of the 
Morse potential, this equation reads 

''"d^ ^^Piv(|/) = 0, y^2rz, 

where F = 2(r — N) + 1. It means that pj^ are the well-known Laguerre polynomials 

(V) 

L}^ {2tz). The following properties of their roots will be used: 

(i) For F > — 1, all roots of L are real positive numbers. 

(ii) As F = —N -I- £, £ — > 0, all the roots of L approach with the asymptotic behavior 



Zj ~ const|£|'''/'^exp (~^^ ■ 

The rational Calogero spin chain with inverse square exchange [27] is obtained as a limit 
of T — > oo, Zj = 1 + T~^/^(^j. The lattice points in this limit are the roots of the Hermite 
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potential Hn{C,)- As F — iV, the lattice becomes equidistant in angles and the model upon 
rescaling is just the trigonometric Haldane-Shastry model [26]. Hence the inhomogeneous 
model defined by the lattice (78a) can be considered as interpolating between Haldane- 
Shastry and Polychronakos model. 

If one chooses as {Pjk} in (77) the spin representation of the permutation group, 
Pjk — (1 + (?j(?A:)/2, the eigenvectors can be treated as in Sections 3-4. Namely, one can 
start from the ferromagnetic vacuum |0 > with all spins up and consider the states with 
given number of down spins M, 

N M 

|V'(^)>= i^{nu-nM)'[l<J:\Q>. 

With the use of the properties of the Lagucrrc polynomials, one finds that in one-magnon 
sector the wave functions can be represented as 



T (r+2m) 

-^iV-l(2T^n) 



The corresponding energies up to universal constant Cn — N{N — l)(3r -|- 2N — l)/24 
are given by 

^i'^ = em = -^(r + m). 

The two-magnon wave functions can be found analytically and the complete set of N{N — 
l)/2 eigenvalues can be written as 

Ejn!n = em + £«(! " 5m,n-i), 0<m<n<A^-l. 
In the M-magnon sector one can find analytically only some eigenstates within the Ansatz 

ip{ni,...nM) = — —F{zn,,...Znj^), 



where F is some symmetric polynomial in {z}. It comprises {N — M + 1)![M!(A'" — 2M -|- 
1)!]~^ eigenvalues which are still additive, 

M 



E{ml} ^^^ruk, TUk+l - 1, < TTlk < N - 1. 



k=l 

This formula allows one to make the hypothesis about structure of the whole set of 
eigenvalues which are described by 

N-l 

Eh. ..Ik — X] ^kh+ii^ — h), 

k=l 

where = —k{r + k)/2 and {Ik} = 0,1. As a consequence of this hypothesis, the 
Hamiltonian H = 2 Yl,^<k hjk(^j(^k is unitary equivalent to the Hamiltonian of the classical 
one-dimensional Ising model with non-uniform magnetic field, 

Ar-2 

(79) 

fe=0 
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with {(Tk} = ±1- This result comprises two above analytical formulae for the spectrum as 
well as Haldane-Shastry and harmonic limits and is confirmed by numerical diagonaliza- 
tion of small lattices up to = 12 with several values of the parameter r. 

The simplicity of the spectrum (79) allows one to compute the free energy / as a 
function of the inverse temperature P in the thermodynamic limit upon rescaling the 
magnon energies with a factor N"'^ [26]. With the use of quasiparticle dispersion law 
e{x) = — x(7 + x)/2 where 7 = T/N one obtains 

f^~^{l ^ d^^og[l + exp(/3e(x))] + j'^ rfa;log[l + exp(-/3e(x))]) , 

which gives at 7 = — 1 the result exactly coinciding with the free energy of the trigono- 
metric Haldane-Shastry model. 

Coming back to the general potential of an external field (76), one has to start with 
the equilibrium equations 

- E "i^'^f + \ E[«^(-. - .7^) + /3 + 7 - (/3 - 7)^7'] = 0- m 

As in the case of the Morse potential described above, let us introduce the polynomial 

N 

Pn{z) = l[{z-Zj) 

with the use of the solutions to (78b) and try to identify the differential equation to which 
this polynomial might satisfy. To do this, note that the function Fj{z)= z{z + Zj){z — 
z j)~^d log p N (z) /dz has simple poles at z = Zk with proper residues, and the equilibrium 
equations can be recast in the form 

TesFj(z)\z=zj = 2aij + Zj{Aa2j - 3a^^) -\- z^(a3j + a^j - 2aija2j) (78c) 

- a^Zj - zj^) + /3 + 7 - (/3 - ^)zf, 

where axj — \p'j^{zj]'~^{d/dz)^'^^pN{z)\z=zy If one supposes that Pn{z) obeys the second- 
order differential equation 

z^p'^ {z)+wi {z)p'^ {z)+W2 {z)pn {z) = (80) 

with some polynomials ^1^2 (-2), one finds upon consecutive differentiations of (80) with 
the use of the formula PN{zj) — that the equilibrium equations in the form (78c) are 
equivalent to 

j-JyW2 + ^-{{a\z' + z-') - z-'wj) + ^-w[ + 2(z((3 + 7) + (/3 - t)^"')] = 0. 
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This condition is satisfied by w;i (2;) = —a{z'^ — l) + {Aa^^P~ji)z,W2{z) = {a — 4f3)z + eN, 
where 71 = 4q;~^7 and parameter e^r is still unknown. One of the solutions to (80) is a 
polynomial of the degree N if the parameters a and (3 are restricted by 

N-1 
(5^ —a. 

The equation (80) is now written as 

z'^p"n{z) - [a{z^ - 1) + (71 + AT - l)z]p'j^{z) + {aNz + e^) = 0. (80a) 

The substitution Pn{z) — z^ -\- Y^^^ ^i^^ results in the recurrence relation for d- coeffi- 
cients in the form 

adi-i{N -/+!) + di[eN + /(/ - 71 - N)] + a{l + l)di+i =0, / = 0, ..N. 

It should be solved under the boundary conditions 

di = 0, (ijv = 1, d^+i — 0- 

The last condition results in A'^th order equation for the parameter e^v- The solution must 
be chosen so as to have all the roots of Pn{z) positive. It is unique since the system of 
particles which repel each other has only one equlibrium point being confined in the field 
with potential (76). 

Due to (80a), various symmetric combinations of the roots of (78b) can be expressed 
analytically in terms of a, 71 and e^. In particular, the energy of classical equilibrium 
configuration does not depend on e^v and is given by 

TV f - 1 9 o\ 

As for corresponding spin chain with the Hamiltonian H = J2j<k hjk{(^j<^k^^)j the strategy 
for finding eigenvalues is the same as for the Morse potential described above. However, 
the information which could be obtained by this way is much more scarce. In M-magnon 
sectors with M < N/2, one can use the ansatz 



V'(ni, ...hm) = — =rp — ri '^Qi^ni, --Z; 



for multimagnon wave function, and show that the eigenequation can be cast in the form 



^(2^' 1,^, d ^ Zjd/dz^-4d/dz, 



M \ 

+ M[(M - 1)(4M + l)/3 - M(7i + A^ - 1) + cat] + a{N - 2M) Y^^k - '^Em\ Q ^ 

k=l J 
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For even N, the solution at M = N/2 [Sz = 0) is given by Q=constant and the corre- 
sponding eigenenergy reads 

En/2 = 1/2{M[(M - 1)(4M + l)/3 - M(7i + AT - 1) + e^]}. 

It was verified numerically that for small lattices {N < 8) at various sets of parameters a 
and 7i this is the exact ground state of the antiferromagnetic chain (77). Unfortunately, 
this approach does not allow to identify other states and write down such a simple formula 
for the whole spectrum as in the case of the Morse potential. 

6. The related Hubbard chains: are they integrable? 

There are another many-body systems on a lattice connected to the Heisenberg-van Vleck 
spin chains discussed above: the itinerant fermions of spin 1/2 which interact being at the 
same lattice site. The corresponding models are Hubbard chains with the Hamiltonian 

HhuO = E tjkct^^ka + 2C/E(4c,-t - l/2)(c+9i - 1/2), (81) 
j¥=k,<7 j 

where the operators c+. create fermion with spin a on the site j, 

{c^a: Cka'} = Sjk5aa': {Cja: Cka'} = 0, (82) 

tjk = t{j—k) is the hopping matrix comprising probability amplitudes for hopping between 
sites j and k (it is supposed to be Hermitian) and C/ > is the strength of on-site repulsion. 

This model was originally introduced by J.Hubbard [42] in three dimensions to describe 
a metal-insulator transition for systems of fermions with spin. It was found that ID version 
(81) is solvable by the Bethe Ansatz [43] in the case of nearest-neighbor hopping under 
periodic boundary conditions, 

tU) = + ^\j\,N-i- (83) 

. The proof of integrability of (81) with the hopping (83), i.e. constructing of the nontrivial 
integrals of motion which commute with (81), came much later [44]. There are two trivial 
invariants: total number of fermions M and number of fermions of up (down) spins which 
are conserved due to su{2) invariance of (81). 

The connection with Heisenberg-van Vleck chains discussed above comes in the limit 
of infinite U a,t M — N (half- filled band). In this limit, fermions are not allowed to 
occupy the site twice and hop, i.e. they can interact only via spin exchange. The spin 
Hamiltonian, which arises in the lowest order in t/U, has the form 

N 

Hspin = E \^jk\ ^j^k- (84) 
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It is this relation on which Gebhard and Ruckenstein (GR) [45] proposed the solvable 
model with hopping 

They were able to guess the simple effective Hamiltonian which comprises all the spectrum 
of Hnuh with hopping (85) but failed in proving this result analytically. Note that till 
now this proof is lacking despite the physical consequences of the GR hypothesis were 
investigated thoroughly [46] and numerical calculations also support it. Moreover, on the 
base of (84) yet another model has been proposed [47] with short-range hopping on the 
infinite lattice, 

t{j) — — isinhre/ sinh(«;j). (86) 

The authors of [47] used the hypothesis of the asymptoic Bethe ansatz for the model 
(86) without any proof of integrability and found quite satisfactory properties in the 
thermodynamic limit. They showed also that (86) includes, as a limit of k ^ oo, the 
nearest-neighbor hopping (83) on the infinite lattice. 

On the base of correspondence with Hnuh ^-nd its limit (84) one can guess also the 
integrability of elliptic model with hopping being some "square root" of elliptic exchange 
(12). But in all these cases, one has to find conserved quantities so as to prove integra- 
bility without appeal to any limit or numerical calculations. This problem is not solved 
completely till now. But some explicit indications to the integrability are found and will 
be discussed later. 

In the spectrum of the model with long-range hopping (85), some degeneracies were 
found similar to the degeneracies for the Haldane-Shastry model [48] . This shows that the 
model might have additional symmetry besides usual one. For Haldane-Shastry model, 
it was found that this symmetry is given by infinite vector algebra, the sl{2) Yangian 
discovered before in [49]. It is natural to try to find at first the source of degeneracies 
for the Gebhard-Ruckenstein model (85). Due to explicit si(2)-invariance of the Hubbard 
Hamiltonian, it is useful to introduce, instead of fcrmion c-opcrators, their bilinear spin- 
like combinations extending the concept of spin to different sites. Namely, the product 
of operators c^^Ckr can be arranged as 2 x 2 matrix {Sjk)'^ labeled by spin indices, which 
allows one to define the S'-operators as 

^fk ~ ^^{'^ "'S'jfc), Sj/. = tT{Sjk), S°' = S^j, Sj = Sjj, 

where a a are the Pauh matrices. Note that Sj'/2 and Sj are the spin density and fermion 
density operators. The commutators of these ^'-operators are 

[^jk^ ^bn — ^klSj^ — SmjS^^, (87) 
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There are a lot of other relations between these operators due to their composite nature. 
Some of them can be written down explicitly, 

Sjk^hn + ^Im^fk + ^Ik^j'm + ^jm^lk — ^jk^tm + ^ImSfk + ^IkSfrn + ^jmS^ki 
^'jk^lm + ^jk^hn + ^'jm^lk + ^jm^U: = '^"^ {^jmi'^^lk ~ 'S'Ifc) + S'J^S'i^ , (88) 

These basic relations contain also a whole list of others which appear upon equating all 
possible combinations of site indices. In terms of S'-operators, the Hubbard Hamiltonian 
reads 

HHut = Y.t^kS% + UY.{{S] - ir - 1/2). (81a) 

The operators of total spin J" = l/^J^jS" commute with (81a), their s/2 commuta- 
tion relations are obtained from (87) by summation over lattice sites, [I°',L^] = ie"^'^!'^. 
Consider now the operator 

J'' -It. {{fjk + hjk{S^ + 5° - 2))Sf, + gjks'^^'^SfS^) (89), 

where fjk = f{j — k) etc. and g and h are odd functions. It is possible to show, with the 
use of (87-88), that Hffub commutes with if the following set of functional equations is 
satisfied [50], 

(gjl - gkl)hjk = -hjlhkl, j T^ky^ly^j, 

lU fjk/2ho + gjkhjk = - t Z] ^ji^ki , j ^ k, 

^ I 

T^fjl^kl — fklhjl) = 0, 

I 

tjk — hohjk, 

where ho is a free parameter. It turns out that the only solutions to these equations just 
give the trigonometric (finite N) and hyperbolic (infinite lattice) forms of hopping (85) 
and (86)! In the trigonometric case one finds 

fjk = 0, gjk = ^ cot(7r(j - k)/N), hjk = i sin-^(7r(j - k)/N), 

whereas in the hyperbolic case 

sinh(K)(j — /c) 1 w /• 7\\ 7 ■ • 1-1/ /• 7\\ 

fjk = u smhUii - k)V Sjk^ ^coth{K{j -k)), hjk^ismh {k{j - k)). 
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Note that J" does not depend on U in the trigonometric case. It is natural to ask of which 
symmetery does this new vector operator correspond. It turns out that this symmetry is 
just Yangian Y{sl2) as it can be seen from the commutation relations 



[J'',K^] + [J^,K''] =0, 



(90) 



where 



and 5—-1 in the trigonometric case and 1 in hyperbolic one. The equation (90) is just 
the defining relation for 5/2 Yangian. Note also that for all odd functions t{j) there is a 
canonical transformation 



which leaves the Hamiltonian invariant but transforms the Yangian generators into 
an independent set of generators /'", J'^ of another representation of SI2 Yangian. It turns 
out that these two representations commute and can be combined to a 1^(5/2) © ^(5/2) 
double Yangian. The fact of this commutativity is nontrivial and is of dynamical origin. 
To verify it and (90), one needs the exphcit form of the functions /, g, h in (89). 

The Yangian operator of the nearest-neighbor chain on an infinite lattice found in 
[51] can be obtained as a hmit of the operator (89) as k ^ 00. In the limit of C/ — > 00 
for half-filled band, where number of fermions coincides with the number of lattice sites, 
one can set Sj — 1 and recover in the trigonometric case the Yangian for the Haldane- 
Shastry model [48]. Thus such rather unhke models as Haldane-Shastry chain and the 
infinite Hubbard chain with nearest-neighbor hopping are in fact connected: they could 
be considered as limiting cases of more general model with the hopping given by elliptic 
functions. 

It is worth noting that the presence of the Yangian symmetry does not imply in- 
tegrability. To prove integrability, one has to construct the set of scalar currents with 
number of its elements at least equal to the number of lattice sites. It was proved for the 
Hubbard model with nearest-neighbor hopping by finding its connection to spin ladder 
and with two coupled six- vertex models [44]. These methods definitely do not work for 
the Gebhard-Ruckcnstein model and its hyperbolic counterpart. One has to find another 
method for constucting integrals of motion. 

To provide examples of the conserved currents which might exist for some choice of 
the hopping matrix, consider the ansatz 




(91) 



JV 



(92) 
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which is most general scalar operator bilinear in {S}. By definition, Ajk = A{j — k) etc. 
The condition [Hnub, ^] = with the use of (87-88) can be cast into the form of two 
functional equations 

AtjkiBik - Bji) + {t.iDik - Djitiu) = 0, (93) 
2{tjkEki + tkjEji) + {tjiDki + tkiDji) = 0, (94) 

the definition of A 

^jk = -2-Djfc + {2U)~^[-8tjkBjk + 2tkjEjk - rjk], 

where 

N 

'"'jk — $Z 'tjlDlki 

and several "boundary" equations for t, B and D: 

N 

X] i'tjl^lk — ^jltlk) = 0, 
N 

'^{tjkDkj — Djktkj) — 0, 

N 

'^^{tjk^kj ~ tkjAjk) = 0. 
k^j 

The first functional equation (93) is just the Calogero-Moser functional equation (10) with 
known general analytic solution. The second functional equation (94) always has solutions 
for Ejk if t and D are given by solutions of (93). Each function in these and "boundary" 
equations can be expressed via basic solution to (93), and the role of "boundary" equations 
is to specify the real period of the corresponding Weierstrass functions, which turns out 
to be N. The basic solution reads 

The other functions in (92-93) are expressed as (recall that tjk = t{j — k) etc.) 

t{x) = tQtl^ix), B{x) = -^^l){x)'il){-x), 

D{x) = d[V''(x)-(^^ + Civ(A) + ^)V'(x)], E{x) = ^^[l-/iV(^ + A)V(-^-A)], 

rix) = todilj{x)[-{N - ^)pn{x) + h^{N - 2)r{x) + {t{x) - h,){2xCN{N/2) - NC{x)) + s], 

N-l 

t{x) = Cn{x + A) - Cn{x) - Civ(A), h = Viv(A)/2, s ^ -{N - 2)p^(A) - Yl Pn{1), 

1=1 

where a^, Cn and pN are the Weierstrass elliptic functions determined by the periods 
cui — N,u!2 — itt/h-i \ = ia ia + N/2, v = K,d,h,a,l3 being arbitrary real 
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parameters. At these conditions, the hopping matrix is Hermitean. Besides this general 
solution, there are the degenerate rational, hyperbolic and trigonometric ones, which 
correspond to one or two periods of the Weierstrass functions. In the first two cases, the 
lattice should be infinite. Checking the absence of "boundary" terms is nontrivial task 
with key formula 

[p{y + A) - p(A)] [Cix -y)- Cix + A) + C(y) + C(A)]+ 

[p{x + A) - prnav - ^) - av + a) + c(^) + c(a)] = p\x). 

These formulas for t, B, D, E, r, A define the scalar current (92) for the model with elliptic 
hopping which comprises all the hopping matrices (83), (85, 86) considered above. At A 
being the half-period of the Weierstrass function, the function iIj{x) becomes odd and 
yet another current is obtained from (92) by the canonical transformation (91). 

The presence of scalar currents commuting with Hamiltonian is the first evidence of 
the integrability for the Hubbard models with the hopping (95) presenting the "square 
root" for the elliptic exchange in Heisenberg-van Vleck chains. It is possible to find the 
corresponding two-fermion function analytically [52]. However, it is seen also that the 
construction of higher scalar currents is extremely hard problem and many-fermion wave 
functions should be also cumbersome and complicated. Till now, nothing is known even 
about ground-state wave function of the simplest trigonometric Gebhard-Ruckenstein 
model: it is neither of Jastrow-typc as for the Haldane-Shastry model nor of Bethe ansatz 
form as in the case of the hopping (83). 

7. Concluding remcirks. 

The main known facts about the integrable Heisenberg-Van Vleck chains with variable 
range exchange and related Hubbard models were reviewed. Many questions in their 
theory are still open. 

As concerns the integrability of these models, understanding it from the Yang-Baxter 
viewpoint is highly desirable. For the spin chains, it is quite probably that the corre- 
sponding R matrix is the same as in [33]. The problem of mutual commutativity of the 
set of operators (18) might be solved in this way. Nothing is known for the integrabihty 
of elliptic Hubbard chains except of the simplest conserved current (92) and two-fermion 
wave function. 

The model with hyperbolic exchange on infinite line should have rich variety of mul- 
timagnon bound states which are given by solutions to transcendental equations 1 — 
{2K)~^[f{pj) — f{pk)] = as it follows from (55). It would be of interest to find more 
simple way of constructing eigenfunctions of the Calogero-Moser Hamiltonian with inverse 
square hypebolic particle interaction. 
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The exact equations of Bethe ansatz type for the case of periodic boundary conditions 
are too complicated at the present stage of finding solutions to quantum elliptic Calogero- 
Moser equation at Z = 1. One cannot exclude the possibility of discovering their more 
simple form which would be of use to verify the hypothesis of asymptotic Bethe ansatz in 
the thermodynamic limit. The construction described above does not allow neither to do 
that nor to establish the correspondence with the trigonometric Haldane-Shastry model. 

In the models on inhomogeneous lattices, the main problem also consists in finding the 
proof of integrabihty for the most general potential of the external field (76). The simple 
formula for the spectrum for the case of the Morse potential, which comprises rational 
and Haldane-Shastry models, still waits analytical confirmation. If one would find the 
explicit form of the unitary transformation of the basic Hamiltonian to its simple effective 
form (79), a lot of results about various correlation functions would be obtained for the 
Haldane-Shastry chain. 

The only known results about the spectra of the related Hubbard models are given 
by original work of Gebhard and Ruckenstein [45]. The trigonometric and hyperbolic 
versions both have the 5/2 ® SI2 Yangian symmetry and scalar integrals of motion (92). 
The most challenging problem is to prove the integrabihty and find the Bethe- ansatz- like 
formulas for the spectrum of the most general Hubbard model with elliptic hopping (95). 
Its solution could clarify the algebraic nature of the integrabihty of all the models under 
present discussion. 
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